Fast domain wall propagation in uniaxial nanowires with transverse fields 
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Under a magnetic field along its axis, domain wall motion in a uniaxial nanowire is much slower 
than in the fully anisotropic case, typically by several orders of magnitude (the square of the di- 
mensionless Gilbert damping parameter) . However, with the addition of a magnetic field transverse 
to the wire, this behaviour is dramatically reversed; up to a critical field strength, analogous to 
the Walker breakdown field, domain walls in uniaxial wires actually propagate faster than in fully 
anisotropic wires. Our conclusions are based on leading-order analytic calculations of the velocity 
and critical field as well as numerical solutions of the Landau-Lipschitz-Gilbert equation. 

PACS numbers: 75.75.-c, 75.78.Fg 



^ ■ The dynamics of magnetic domain walls in ferromagnetic nanowires under external magnetic fields [lHlO( and 

spin-polarised currents [Iol - [l5j is a central problem in micromagnetics and spintronics, both as a basic physical 
phenomenon as well as a cornerstone of magnetic memory and logic technology @, [l3|-[l5[ • From the point of 
view of applications, it is desirable to maximise the domain wall velocity in order to optimise switching and 
response times. 

Partly because of fabrication techniques, attention has been focused on nanowires with large cross-sectional 
aspect ratio, typically of rectangular cross-section. In this case, even if the bulk material is isotropic (e.g., 
permalloy), the domain geometry induces a fully anisotropic magnetic permeability tensor, with easy axis along 
the wire and hard axis along its shortest dimension (lol ;23j . Nanowires with uniaxial permeability, characteristic 
of more symmetrical cross-sectional geometries (e.g., square or circular), have been less studied [l7l - [l9j |. Here 
we report that, in the presence of transverse fields, the domain wall velocity in uniaxial wires is actually faster 
than in fully anisotropic wires up to a certain critical longitudinal field. This is in sharp contrast to the case of 
zero transverse field, where the uniaxial velocity is much less than in the anisotropic case. 

We employ a continuum description of the magnetisation. For a thin nanowire, this is provided by the 
one-dimensional Landau-Lifshitz-Gilbert (LLG) equation (2p| - [23l |. which we write in the non-dimensionalised 
<— | . form 

q . M = M x H - aM x (M x H) . (1) 

Here M.(x,t) is a unit-vector field specifying the orientation of the magnetisation, which we shall also write in 
polar form M = (cos 0, sin 6 cos $, sin 6 sin $). The effective magnetic field, H(m), is given by 

> ■ 

ON ■ H = Am" + Kmux - K 2 m 2 y + H a , (2) 

where A = 1 is the exchange constant, K\ = 1 is the easy-axis anisotropy, K 2 > is the hard- axis anisotropy, 
. H a is the applied magnetic field (taken to be constant), and a is the Gilbert damping parameter. Domains 
' correspond to locally uniform configurations in which M is aligned along one of the local minima, denoted m + 
, and m_, of the potential energy 

[/(m) = -i(m 2 -^ 2 m 2 )-m-H a . (3) 



Two distinct domains separated by a domain wall are described by the boundary conditions M(±oo, t) = m±. 

For purely longitudinal fields H a = Hik and for Hi below the Walker breakdown field Hyy = aK 2 /2, 
the domain wall propagates as a travelling wave [![, the so-called Walker solution <d(x,t) = 6w{x — Vwt), 
d , $(x,t) — 4>w, where 9w and 4>w ar e given by 

9 W (0 = 2 tan" 1 (e" 5/7 ) , sin 20^ = H l /H w . (4) 

The width of the domain wall, 7, is given by 7 = (1 + K 2 cos 2 (frw)' 1 ' 2 , and the velocity is given by 

V w = -l(a + l/a)H 1 . (5) 

For Hi > flwjthe domain wall undergoes non-uniform precession and translation, with mean velocity decreasing 
with Hi P, The effects of additional transverse fields have been examined recently [1, Q. 

If the cross-sectional geometry is sufficiently symmetrical (e.g., square or circular), the permeability tensor 
becomes uniaxial, so that K 2 = [16l. |23||. The dynamics in this case is strikingly different. The LLG equation 
has an exact solution, &(x, t) = 9q{x — Vpt), <&(x, t) = —Hit, in which the domain wall propagates with velocity 

V P = -aHi (6) 
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and precesses about the easy axis with angular velocity —Hi [IT], [3 ■ The precessing solution persists for all 
Hi - there is no breakdown field - but becomes unstable for Hi > 1/2 [l9j . 

For Hi < Hw, the ratio Vw/Vp is given by j{a~ 2 + 1). For typical values of a (0.01 - 0.1), the uniaxial 
velocity Vp is less than the fully anisotropic velocity Vw by several orders of magnitude. As we show below, 
applying a transverse field Hi > to a uniaxial wire dramatically changes its response to an applied longitudinal 
field H\. The transverse field, analogous to hard-axis anisotropy, inhibits precession and facilitates fast domain 
wall propagation. For Hi less than a critical field Hi c (H 2 ), given by ([T8jl below, which depends on H 2 , there 
appears a uniaxial travelling wave in which the domain wall velocity exceeds that of anisotropic wires. 

We first obtain a general identity, of independent interest, which relates the velocity of a travelling wave 
M(ac, i) — m(x — Vt) (assuming one exists) to the change in potential energy across the profile (for zero 
transverse field, this coincides with results of Q and 0). Noting that M = — Via.', we take the square of ((T|) 
and integrate over the length of the wire to obtain 

y 2 ||m'|| 2 = (l + a 2 )||m x H|| 2 . (7) 

Here we use the notation 

/oo 
U'vdi (8) 
-oo 

for the L 2 -norm and inner product of vector fields (analogous notation for scalar fields is used below). Next, 
we take the inner product of Q} with H to obtain 

V(m',H.) = -a\\m x H|| 2 . (9) 
Noting that m' • H = (im' ■ m' — E7(m)) , we combine ([7]) and © to obtain 

V = \{a + l/a) ||m'||- 2 (f/(m_) - C/(m+)) . (10) 

The identity (|10p has a simple physical interpretation; the velocity is proportional to the potential energy 
difference across the wire, and inversely proportional to the exchange energy of the profile. 

From now on, we consider the uniaxial case K 2 = and applied field with longitudinal and transverse 
components Hi,H 2 > (by symmetry, we can assume H3 = 0) with |H a | < 1. An immediate consequence 
of (|10p is that, in the uniaxial case, the velocity must vanish as Hi goes to zero. For when Hi = 0, the local 
minima m± are related by reflection through the 23-plane, and U(m+) = J7(m_). 

To proceed, we substitute the travelling wave ansatz @(x, t) = 9{x — Vt), &(x,t) — <f){x — Vt) into the LLG 
equation fl} to obtain the system 

0" - | sin 2(9(1 + (t>' 2 )- Hi sin 9 + H 2 cos 9 cos 

= -V(a9' +sin00')/(l + a 2 ), (11a) 
sin 6 4>" + 2 cos 6 O'tf - H 2 sin <j> 

= V{9' -asmd<P')/{l + a 2 ). (lib) 

The boundary conditions 9±, 4>± at £ — ±cxd are given by 

sin6» ± cos0 ± + His\n9 ± - H 2 cos9± = 0, cf> ± = 0. (12) 

The system (|11[) is nonlinear, and exact solutions are not known (except, for H 2 = 0, the precessing solution). 
A widely used scheme, due to Slonczewski (23 . [25l |. gives approximate expressions for the position and orientation 
of the centre of the domain wall. Here we use a different approach, valid when either Hi or H 2 are small. The 
simplest case, which we consider first, is where both are small. We write Hi — ehi, H 2 = eh 2 , V = ev, and 
expand the profile to first order in e as 9 = 9q + e9i , <j> = <fio + e</>i ■ Substituting in (jlip we obtain at order e° 
the equations 

9' Q = — sin 9q, <po — constant, (13) 

which characterise the family of zero-field uniaxial optimal profiles with boundary conditions do- = ft and 
$o+ = 0, with solution tan(# (£)/2) = exp(— £). 

At order e, it is convenient to introduce ui = sin^o^i: which, in contrast to 4>i, is well behaved as 9q 0. 
We obtain the following linear inhomogeneous equations for 9i and Ui m . 

L9i = (a + l/a)v9' — hi sin^o + h 2 cos^o cos^o, (14a) 
Lui — — (a + l/a)v9o — h 2 sin^o, (14b) 
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where 

j2 qui 

L = -^2+^ W^) = ^r = l-2sech 2 £ (15) 

and 6\, u\ are required to vanish as £ — > ±00. 

The velocity v and orientation </>o are determined by the solvability conditions for (|14[) . According to the 
Fredholm alternative, given a self-adjoint operator L on L 2 (R), a necessary condition for the equation Lf = g 
to have a solution / is that g be orthogonal to the kernel of L. If this is the case, a sufficient condition is that 
the spectrum of L is isolated away from 0. (In fact, some additional care is required here, as g does not vanish 
as £ — > ±00, but this poses no new essential difficulties.) From (fl~5|) it is clear that 9' belongs to the kernel of L, 
and since the eigenvalues of a one-dimensional Schrodinger operator are nondegenerate, it follows that 9' spans 
the kernel of L. Since W(£) — > 1 as £ — > ±00, it follows that spectrum of L is discrete near 0. (In fact, W(£) is 
a special case of the exactly solvable Poschl- Teller potential, but we won't use this explicitly.) 

Requiring the right-hand side of (|14ap to be orthogonal to 9' and noting that (0' ,6' ) = 2, (9' ,sin6o) = —2, 
and (9' 0i cos9q) = yields 

V = -(a+ l/a)H x . (16) 

Thus, the velocity is independent of H2 to leading order (but see (|2"3"|) below). It is readily checked that (TlG]) 
is consistent with the identity (fTO)) . Requiring that the right-hand side of (|14b|) be orthogonal to 9' and using 
(|TB)) . we obtain 



Requiring sin^o < 1, we obtain the critical field 

H lc = a^H 2 . (18) 

We can extend the calculation of the velocity to the regime where H2 is no longer small. We continue to take 
Hi = ehi and V = ev and expand the profile to first order in e as above. The e°-equations become 

9' = (H 2 -sm9 ), O =O, (19) 

with boundary conditions sin#o± = H2, #0+ > tt/2 and #0- < 7r/2. Thus, for finite H2, azimuthal symmetry 
is broken at lowest order, and the static profile aligns along the transverse field (the alternative, </>o = 7r, also 
satisfies the e°-equations but has higher energy and is unstable to perturbations). The static profile is given by 



"0 K 
tan — = tanh 

2 m 



h (20) 



where k = yl - H\ . At order e we obtain the linear inhomogeneous equations 



a 



L8 X = - 7;v9'r, - hismOo, (21a) 

l + a 2 u ^ ' 

Mq> 1 = —}—v(cos9o)' 1 (21b) 



where 



j2 qui j j 

L —^ + t' M = -^ Shl2 ^ + ff2Sin ^ (22) 

9q is given by (|20[) . and 9±, <f>\ are required to vanish as £ — > ±00. As above, we may apply the Fredholm 
alternative to determine the solvability conditions for (|21[) . One can show that the kernel of L is spanned by 
9' Q and that its spectrum is bounded away from zero. Requiring 9' to be orthogonal to the right-hand side of 
(|21a[) , we obtain 

V = -(a+ l/a)(l - (H 2 /k) cos" 1 H 2 y 1 H 1 . (23) 

It is readily checked that (|2"3"]l coincides with (Ti~6"|) for H2 = and therefore gives ^-nonlinear corrections to 
the velocity, and moreover is consistent with the general identity (jlOp . One can also show that M has trivial 
kernel with spectrum bounded away from zero, so that (I21b|) is automatically solvable. 
From © and 

VW/y = 7 (l-(i7 2 / K )cos- 1 J ff 2 ) < 1. (24) 
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Thus, to leading order in Hi, the domain wall velocity in a uniaxial wire with transverse field exceeds the 
Walker velocity. Numerical results below establish that this continues to hold up as Hi approaches the critical 
field H u . 

To verify our analytical results, we solve the LLG equation ([T]) using a finite-difference scheme on a domain 
—L < x < L where L = 100 (the domain wall has width of order 1). Neumann boundary conditions, m' = 0, 
are maintained at the endpoints. The damping parameter a is taken to be 0.1 throughout. As initial condition 
we take the stationary profile, with 8q given by (|20p and <f>Q — 0. After an initial transient period, during 
which the asymptotic values of m at x — > ±L converge to m±, a stable solution emerges, in which the domain 
wall propagates with a characteristic mean velocity V. (For convenience, we have taken Hi < 0, so that V is 
positive.) In Figure [TJ V is plotted as a function of |i?i| for three fixed values of the transverse field: H 2 = 0.2, 
H 2 =0.1, and the limiting case Hi = 0, where the dynamics is given by the precessing solution. For nonzero 
H 2l the velocity exhibits a peak at a critical field \Hi c \, which depends on H 2 . Figure [2] shows the dependence 
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FIG. 1: Average propagation velocity I 7 of a domain wall as a function of the driving field \Hi\ for three values of 
the transverse field Hi- Inset: The numerically computed velocity (solid curves) and the analytical expression Eq. (123[) 
(dashed lines). 

of |i/i,c| on H 2 , in close agreement with the analytic result ((T8j) . 
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FIG. 2: The critical driving field |i2i, | as a function of the transverse field H2. A linear fit (blue curve) through the 
numerically computed data (blue diamonds) is plotted alongside the analytical result (|18p (red curve). 

As in the fully anisotropic case, the properties of the propagating solution are qualitatively different for driving 
fields \H\\ below and above the critical field. This is demonstrated in Figure [31 which shows contour plots of 
the magnetization in the (x, i)-plane. Figs.|3Ja) and[3Jb), where Hi = —0.01, exemplify the case \Hi\ < \Hi c \. 
The magnetisation evolves as a fixed profile translating rigidly with velocity V. For \H±\ > \Hi c \, as exemplified 
by Figs. [3jc) and[3jd), where Hi = —0.05. the magnetization profile exhibits a non-uniform precession as it 
propagates along the nanowire, with mean velocity decreasing with H±. For \Hi\ 3> Hi Cl the solution can be 
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regarded as a perturbation of the precessing solution. 

In summary, we have established, both analytically in leading-order asymptotics and numerically, the existence 
of travelling wave solutions of the LLG equation in uniaxial wires in applied fields with longitudinal and 
transverse components. We have obtained analytic expressions for the velocity, (|16D and (j2"3")l . and for the critical 
longitudinal field, (fT8|) . above which the travelling wave solution ceases to exist. The analytic expressions are 
confirmed by the numerics. Our principle result, (pHf. establishes that below the critical field, the domain wall 
velocity in uniaxial wires exceeds that in fully anisotropic wires. It would be desirable to extend the asymptotic 
analysis developed here to general anisotropy and to applied fields which exceed the critical field. We hope to 
address these problems in future. 
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